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Abstract. We present a detailed study of the ionisation probability of H and HJ 
induced by an short intense laser pulse. Starting from a Coulomb- Volkov description 
of the process we derive a multipole-like expansion where the first order transition 
matrix may be factored into two contributions: one that accounts for the effect of the 
C^ ' electromagnetic field on the free-electron final-state and a second factor that depends 

lyj , only on the target structure. Such a separation may be valuable to solve complex 

atomic or molecular systems as well as to interpret the dynamics of the process in 



£^ ■ simpler terms. We show that the series expansion converges rapidly, and thus the 

(— I I inclusion of the first few terms are sufficient to produce accurate results. 
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1. Introduction 

The interaction of strong and short laser pulses with atoms and molecules has received 
renewed attention in the recent past, both experimentally and theoretically [1-5], mainly 
because the advances in laser technology have made feasible time-resolved measurements 
of atomic and molecular processes. These advances made possible new experimental 
investigations of atomic and molecular processes on an ultrashort time-scale, and under 
ultra intense laser radiation. These techniques lead the way to measurement of highly 
nonlinear phenomena, and even their control is now possible [6-8]. 

On the theoretical side, since the first works by Kulander [9, 10], more than two 
decades ago, many computational techniques have been developed to solve the three- 
dimensional time-dependent Schrodinger equation (TDSE) for single-electron systems 
[11-14]. The advances in computing power allows us to perform these computations in 
a question of minutes nowadays. However, the computation power demands imposed 
by precision numerical solutions of the Schrodinger equation quickly scales out-of-reach 
when the complexity of the systems increases. 

Alternative, complementary approaches based on perturbative expansions could 
provide results at relatively low-computer cost at the expense of some precision-loss. 
Traditionally, the Keldysh-Faisal-Reiss or Strong-Field Approximation (SPA) [15-17] 
(see also [3]) has been employed to theoretically describe these processes when the field 
is intense. On the contrary, when the field is weak, a First Order Perturbation theory 
(FPA) is sufficient to describe the process [18]. In the intermediate regime, where the 
interaction of the electrons with the nucleus and with the electromagnetic (EM) field 
are of similar strength, none of the potentials may be neglected. To describe ionisation 
processes in this regime was developed the Coulomb- Volkov approximation (frequently 
called CV2 or CVA) [19-22], where both interactions with the ejected electron are taken 
into account at the same level in the final-state. 

The CVA approach has been widely used to investigate the ionisation of atomic 
hydrogen [22-24], alkali metal atoms [25], simple molecules [2G] and positronium [27]. 
Furthermore, there are several implementations of the Coulomb Volkov approximation, 
that may be used in a wide range of conditions: the renormalised CV (RCV) extends 
the CVA to the non-perturbative regime [28], the modified CV (MCV) introduces the 
coupling to intermediate bound states [29], and the doubly distorted CV (DDCV) [30, 31] 
includes the distortion by the laser field also in the initial state. 

Alternatively, Dimitrovski et. al. [32, 33] have presented analytical formulas for 
ionisation by very short pulses, that are independent of the laser intensity, in the context 
of First Magnus Approximation (FMA). When weak and short field are applied, FPA 
and FMA can be used to obtain the sudden approximation [18, 32] where the ionisation 
amplitude is proportional to the momentum transfer and the dipole transition matrix. 

The interaction with intense laser fields could induce, besides multiphoton 
ionisation (MPI), another important phenomenon: High order Harmonic Generation 
(HHG). In a semiclassical recoUision picture it proceeds in three steps: ffist the electron 
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is released from the target, then the electron interacts with the laser field, and in a 
third step, if the electron returns to the atom, radiative recombination may take place. 
Alternatively, in the third step, an elastic scattering or an (e,2e) process may take 
place, giving place to High-order Above- Threshold Ionisation (HATI) or Non-Sequential 
Double Ionisation (NSDI), respectively [5, 34]. 

Despite the well-known shortcomings of the Strong Field Approximation (SEA), 
such as the discrepancies for calculations on different gauges, three-steps processes may 
be understood within its simple framework. The first-order SFA is able to describe 
'direct' ionisation, but a second order term is necessary to include the post-ionisation 
interaction between the free electron and the target. In contrast, this mechanism 
is already included in the first-order CVA, providing both direct and rescattering 
amplitudes [3]. Recently, an alternative approach to this problem, called Quantitative 
Rescattering theory (QRS), was developed by Lin et. al. [34]. In this approach the yields 
for HHG, HATI and NSDI can be expressed as the product of the returning electron 
wave packet probability with photo-recombination, elastic electron scattering, and 
electron-impact ionisation cross sections, respectively. Also recently, in a full quantum 
mechanical description, Frolov et. al. have deduced an analytical factorisation of the 
spectra in terms of an electron wavepacket and the cross-section of photo-recombination 
(for HHG), and elastic electron scattering (for HATI), see [35-37] and references therein. 

This type of yield factorisation allows the study of more complex atoms or 
molecules, by separating the roles of the laser-pulse and the parent-ion, and allowing 
to accurately extract information about the target, (see for example [34, 38-40]) 
offering a promising tool for dynamic chemical imaging with temporal resolution of 
few femtoseconds. 

In this work we study the above threshold ionisation (ATI) of atomic hydrogen by 
short laser pulses within the CVA approach. We derive an expansion of the transition 
matrix in powers of the laser-field vector potential, whose first term is sufficient to 
describe ionisation in the multiphotonic regime. This first-order term is written as the 
product of an one-photon ionisation transition-matrix and an integral factor depending 
on the laser pulse. 

We compare the ATI spectra with TDSE results and analyse the convergence of the 
expansion, how each factor of this approximation contributes to the electronic spectra, 
and present some useful applications of this factorisation. Finally, we also investigate 
the ionisation from some excited states of the hydrogen atom, and from Hj employing 
exact wavefunctions for both the initial and final channels. 

Atomic units are employed, except where otherwise stated. 

2. Review of the Coulomb- Volkov approximation 

Let us consider the ionisation of an atomic or molecular system by interaction with a 
finite laser pulse of duration r. In the length gauge, the action of the pulse may be 
described as a time-dependent force produced by the electric field. Thus, the time- 
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dependent Schrodinger equation (TDSE) reads 



.d^{r,t) 



P^ 



+ V{r,t) + F{t)-r 



^(r,t) 



(1) 



We consider a finite laser pulse, that exerts a force on the system given by its 
electric field 



F{t) = Fo sin {uj{t - to)) sin^ (7rt/r) i 



(2) 



for < t < r, and that vanishes at all other times. The pulse is characterized by its 
central frequency u, the pulse duration r, the polarisation vector i, and the phase-shift, 
chosen as wto = {^t — 7r)/2 for symmetric pulses. 

While the perturbation vanishes outside the time-interval AT = (0, r) we can write 
the ionisation transition-matrix, in its prior form, as 



T, 



fi 



(vI/7(t)|F(t)-r|0,(t))dt. 



(3) 



Here "^jit) is the exact wavefunction for the final state with ingoing boundary 
conditions, and 0j(t) is the asymptotic target wavefunction in absence of external fields. 
The double differential ionisation probability in energy and angle of the emitted 
electron is obtained from the transition matrix magnitude as 

The Coulomb- Volkov approximation (CVA) is obtained by replacing the exact final 
wavefunction 1^7(7", t)) by a product of factors corresponding to the solution of two 
separated problems: one for the isolated field-free atom and one containing the effects 
of the electron evolving in the external electromagnetic (EM) field. In order to get 
the correct asymptotic behaviour the plane- wave part must be corrected [19, 20]. This 
procedure is similar to the one carried-out in continuum-distorted-wave (CDW) theories 
developed for ion-atom collisions many years ago [41-43]. In this approximation the 
final-state wavefunction takes the form [19, 20] 

„ifcr 
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A(t') dt' 



{A{t')f dt' 



A{t) 



Fit') dt' 



(5 a) 
(56) 

(5c) 
ihd) 



Here y? r (r) is the final-state of the electron ejected with momentum fc, corresponding to 
energy Ef = k'^/2, and with ingoing wave boundary conditions. We have explicitly set 
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m = 1 and Z = —1 for the electron mass and charge, respectively. For a pure Coulomb 
potential interaction (hydrogen atom) the distortion factor takes the familiar form: 

Dc{k, r) = N~{u) iFi (iz/; 1; -i(/cr + k ■ r)) , 

where the normalization factor is defined is given by N~{u) = r(l — ii/)e~'^'^/^, and 
u = —1/k is the Sommerfeld parameter. For general, more complex, targets the 
distortion factor must be obtained numerically, and is defined in terms of the final- 
state target eigenfunction ipj as given by {5b). 

After some manipulation, the CVA transition matrix may be written as [21, 22] 



T^' = f{0)g{0)+ dth{t)f{t)g{t) (6a) 



where 



h{t) = i (w/, + k ■ A{t) + A'^{t)/2) , {6b) 

f{t) =exp{j\{t')dt'), (6c) 

g{t) = JdripJ*{r)exp{-iA{t)-r)ip,{r), {6d) 

and Ufi = Ef — Ei. For symmetric pulses A{t) = 0, and an additional term involving 
g{T) in the above expression, vanishes due to the orthogonality of the initial and final 
wavefunctions. 

This approximation has been successfully employed to describe atomic ionisation 
by short UV laser pulses. However its validity is limited to the perturbative regime, 
where the probability of ionisation is small. It has been shown that the approximation 
fails when the photon energy is smaller than the ionisation potential u <^ Ip. In order to 
extend the range of applicability of the theory it has been developed a modified version 
that includes intermediate bound states [24, 29]. 

3. Multipolar expansion of the T-matrix 

It has been shown that the probability may be evaluated efficiently, rendering good 
results with relatively low computational cost in the case of atomic hydrogen. However, 
its extension to more complex systems such as many-electron atoms or even small 
molecular ions can considerably increase computation times. It is desirable to develop 
approximations that keep the good performance of the CVA but may be applied to 
laser-induced ionisation of large molecules. 

In the range of applicability of the CVA theory, the amplitude of the field Fq is 
small and u > Ip. Then, in most cases the magnitude of the vector field A{t) is also 
small. In that case we may expand the exponential in g{t) in (6c/) and keep only the 
first, most significant, terms 

exp (-iA ■ r) = 1 - iA ■ r {A ■ r)^ + . . . 
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We obtain a series expansion for the transition matrix T9y = T^^'> +T^^'^ + T^'^'> + . . ., 
whose first terms are given by 

tW = (^7|^,) (7a) 



T(i) = -iL«. 



rpi2) ^ _1^(2) 

2 



(76) 
(7c) 



/(0)A(0)+ / dtf{t)h{t)A{t) 
Jo 

f{0)A\0)+ f dtf{t)h{t)A\t) 
Jo 

where we have defined the matrix elements L^^^ = ((/)7|r|y9j) and L^"^^ = {Lpj\(e ■rY\(pi). 

Keeping only the first non- vanishing order allows us to decouple the results as two 

factors from very different origin, one of them accounts for the laser perturbation and 

its effect on the final state while the other is sensitive to the target structure: 

T(i) = -i X(^) • M(^) , (8) 

where X^^^ is the one-photon ionisation transition matrix, and 

mW = f F{t)f{t)dt. (9) 

Jo 

The above approximation relies on the condition A{t) ■ r ^ 1, for that reason 
we call it dipole approximation or DipA. While the amplitude of the vector potential 
A{t) is proportional to the electric field amplitude Fq and inversely proportional to the 
frequency w, the range of validity of the approximation is constrained by these two 
parameters. In principle some of these constrains are already considered in Coulomb- 
Volkov approximations. It has been shown that the CV2 approximation provides 
reliable results of atomic and molecular ionisation by ultraviolet lasers whenever the 
laser intensity is compatible with the perturbation regime and the photon energy is 
higher than the electron binding energy [22, 29]. However, for a given EM field this 
approximation may fail as highly-excited states are considered, because the radius of 
the atomic or molecular system increases. We note that this approximation may also 
be obtained by neglecting the first factor in the Volkov state (5c). Because this factor 
contains the spatial depencency of the EM field, space and time are then decoupled in 
the transition matrix integrals (3). 

The Keldysh parameter F = ^J{Ip/2Up) ~ 1 defines the limit between 
multiphotonic and tunneling regimes for ionisation processes. Here the ionisation 
potential energy is Ip and the ponderomotive energy Up = F^/Au^. Then we expect 
that DipA will work well in the multiphotonic regime, when F ^ 1, since in this case 
A ~ 1/F will be small. 

To analyse the spectra obtained in DipA we study the different factors in (8). As we 
mentioned before, X^^^ is the ionisation transition matrix due to the absorption of one 
fictitious photon with frequency Ufi = Ef — Ei. For eigenf unctions of the unperturbed 
Hamiltonian, the matrix element can be written equivalently either in the length or 
velocity gauge: 

xW = (^^|r|v9,) = - — (^7|V|(^,) (10) 
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This term does not contain any information on the EM field, and will be the same 
for all laser pulses. On the other hand, since the function / can be factorised as 
/(t) = j{t) exp (iufit), the contribution of the field can be written as 

mW = f dtF{t)j{t)e'^f'' (11) 

Jo 

Here j{t) is the part of the Volkov state that is taken into account in DipA and a phase 

that could be omitted: 

j{t) = e-'^-f'^ exp 



ik- f A{t')dt' 



X exp 



^ / iA{t')f dt' 



(12) 



Equation (11) have the aspect of a Fourier transform. However, the function j 
and the Fourier frequency Ufi are not independent, since both depends on the electron 
energy Ef. For sufficiently low intensity of the field, j ~ 1 and this equation defines 
exactly the Fourier transform of the EM pulse. Also, we note that the only information 
about the target that TVf ^^^ contains is the energy Ei in uji. For this reason, spectra 
for different targets, or initial states with different binding energies can be reproduced 
approximately, by shifting the final electron energy. 

In the following sections we analyse the present multipolar expansion. We present 
the CVA spectra for H atoms and H^ molecular ions, comparing with TDSE results 
and the spectra obtained employing the first- and second-order transition matrix in the 
multipolar expansion of CVA. 

4. Spectra for H 

4.1. Comparison of different theories 

In this section we compare CVA and DipA spectra to the results obtained solving the 
TDSE with the code QPROP [12]. For CVA calculations we employed equations (4) 
and (6 a), while in the DipA case we considered (4) and (8). 

In figures 1, 2 and 3 we present the density of probability (DOS) for the ionisation 
of atomic hydrogen (Is). The spectra, as function of the electron energy, were obtained 
integrating over all electron-emission directions. We consider laser pulses of intensity 
Fq = 0.05 a.u., duration corresponding to N = 1, 7 and 27 cycles, and frequencies u = 
0.855 (Fig. 1), 1.71 (Fig. 2), and 0.427 a.u. (Fig. 3). The results for N = 1 and 27 are 
multiplied by 10 and 10~^, respectively, for better visualisation. 

In general, we observe from these figures a good agreement between CVA, TDSE 
and DipA results, specially for high laser frequencies. These examples are within the 
limits of applicability of the CVA approximation; thus they are almost indistinguishable 
from the numerically exact results given labelled TDSE. At lower frequencies the 
disagreement between CVA and TDSE has been analysed by Duchateau et. al. [21, 22], 
and they pointed out that the CVA fails because it does not consider transitions to 
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Figure 1. (colour online) Electron spectra for ionisation of H(ls) as function of the 
ejected electron energy for a laser pulse with N = 1 (up), N = 7 (middle) and N = 
27 (bottom) cycles. Laser frequency a; = 0.855 a.u., Fq — 0.05 a.u. (F = 17.1). Full 
line: CVA, dotted line: Dipolc Approximation (DipA) and full line with circles: TDSE 
results. 
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Figure 2. (colour online) Idem Fig. 1 with w = 1.71a.u. (F = 34.2) 



intermediate excited states. This mechanism is more important when the ionisation 
energy is larger than the laser frequency. 

Also, we observe the formation of above-threshold-ionisation (ATI) peaks when 
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Figure 3. (colour online) Idem Fig. 1 with a; = 0.427a.u. (F = 8.55) 

increasing the number of cycles of the pulse. These structures correspond to the 
absorption of n photons at the value energy En = Ei + nu — Up (see for example 

[3]). 



4.2. Analysis of Dip A 

Let us now discuss the DipA results, plotted with dotted line in figures 1 to 3. In 
these figures we observe that the shape of the spectra is well reproduced by the dipole 
approximation, but there are disagreements in the ATI peaks, corresponding to multiple- 
photon absorption, where DipA underestimates the ionisation probabilities. In all cases, 
the first peak, that is the most important in magnitude, is well reproduced by the 
DipA. As a consequence, total ionisation probabilities are accurately given by this 
approximation. This is observed in figure 4, where we show total ionisation probabilities 
for a frequency uj = 0.855 a.u, as function of the laser pulse duration, and for several 
EM intensities. We can observe that the first-order approximation (DipA) reproduces 
very well the full calculations in a very extended range of laser pulse parameters. 

As we mentioned before, the intensities of the secondary ATI peaks are not well 
described by DipA. This means that in this regions there are couplings between time and 
space, and/or the exponential factor exp (— iA • r) plays a crucial role in the description 
of the ionisation spectra. 

The u = 0.427 a.u. case of figure 3 is a bit different since the ATI peaks are closer 
to each other and shifted to the left. The first peak would appear at negative energies, 
below the ionisation threshold. Thus, the first peak corresponds to the absorption of 
two photons, showing worse agreement. This is an expected result, according with the 
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Figure 4. Comparison of hydrogen (Is) total ionisation probabilities by pulses laser 
of frequency lo — 0.855 a.u, for several different intensities. Lines correspond to full 
CVA calculations, while symbols were obtained in the first order DipA. 



above description, for low frequencies or high amphtudes A. 

We note that for large electron energies, outside the ATI peaks, the background 
of the spectra is also well reproduced. This fact can be explained noting that for high 
energies the contribution of the term k- j A{t) dt is more important than the contribution 
of r ■ A. Then, this last term can be omitted in the exponential of the Volkov state (5 c), 
resulting precisely in DipA. 

Surprisingly, short pulses corresponding to A^ = 1 cycle are well reproduced by 
DipA independently of the laser frequency (figures 1 to 3). For very short pulses and 
in the context of the First Magnus Approximation (FMA), Dimitrovski et al. [.32, 33] 
found that the ionisation amplitude is independent of the pulse shape and, therefore, 
all pulses are equivalent if they have the same value of A(0). In connection with the 
present Dipole Approximation, we note that pulse amplitudes A(t) that vary slowly in 
time allow naturally the decoupling of spatial and temporal integrals in the transition 
matrix. 

Furthermore, when for short pulses and weak fields, FMA and first order 
perturbation theory (FPA) can be combined to give rise to the sudden approximation 
[18, 32], where the ionisation amplitude is described by (8), but simplifying /(t) = 1 
in (9). Our proposed DipA presents the same type of benefits that the sudden 
approximation, but for a extended range of laser parameters, as we can observe in 
figure 4. 
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4-3. Dip A applications 

DipA enable us to understand the features of the spectra, allowing us to separate those 
that come from the pulse from those that are exclusively due to the nature of the 
target. For example, because we consider ionisation of H atoms from the ground state 
in the cases presented in figures 1 to 4, in all cases the contribution of L^^' is the same. 
Its squared modulus, integrated in electron emission direction, is proportional to the 
photoionisation cross-section and is a monotonically decreasing function of the energy 
(see for example eq. 37 of [11]). 

On the contrary, the factor M^^\ is different for each of the nine spectra plotted 
in figure 1 to 3, corresponding to the nine different laser pulses. For symmetric pulses 
with A^ = 1 cycle, the shape of the pulse is basically a peak, then the Fourier transform 
predicts a widely spread spectra, as observed from the above figures. Increasing the 
number of cycles, the Fourier transform of the pulse narrows considerably, as a delta- 
like function, evaluated in the central frequency of the pulse. This behavior gives rise 
to the first ATI peak. To analyse the formation of the others ATI peaks we can expand 
the exponential in j{t) in powers of the argument. Each of these terms contributes 
to an harmonic of the original pulse, for this reason his Fourier transform shows the 
ATI peaks centered in positions that are multiples of the central frequency u. Also, 
when increasing the intensity of the pulse, more terms in the exponential expansion are 
needed, explaining the formation of new ATI peaks as the intensity increases (see for 
example figure 1 of [22]). 

Summarising, we can note that the general shape of the spectra is defined by M'^^'^: 
the number of ATI peaks, their amplitude and width could be estimated from the 
analysis of this factor, as well as the dependence on the carrier envelope phase (CEP), 
fixed to cjto iii this work. 

One important application of this approximations is its use in more complex 
systems, like molecules. In this case it is possible to construct the spectra for laser pulse 
ionisation by multiplying the one-photon ionisation transition-matrix of this complex 
system by the factor M'^^^ with the information of the pulse. For sufficiently high- 
frequency pulses, as the generated from typical HHG spectra [15-48], this approximation 
will be valid. Conversely, knowing the emission spectra of an unknown target interacting 
with such pulses it is possible to decompose the spectra, as proposed in DipA, and obtain 
information about the target. This type of experiment where the extreme-ultraviolet 
attosecond light-pulse are used as source to emit electrons and obtain "tomographic 
images" of the target have been achieved recently [8, 49, 50]. 

Another interesting application of the DipA is its application to a simpler Strong 
Field Approximation (SEA). The first order of SEA in length gauge (LG SEA) can be 
obtained replacing the coulomb wave function (fl by a plane-wave in the present CVA 
calculations. It is known that this approximation may overestimate the exact results 
and needs to be renormalised [21]. However, it is possible to employ SEA with the 
DipA factorisation in the regime where DipA works well (F ^ 1 and u > Ip). Since 
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Figure 5. (colour online) Ionisation spectra for H(ls) for laser pulses with cj = 0.855 
(up) and 1.71 a.u. (down), N = 7 cycles and Ff, = 0.05 a.u. Full line (red): DipA-GV 
(SFA) see text. Dotted line (blue): first order LG SFA. Full line with circles (black): 
TDSE results. 



the plane-wave is not an eigenfunction of the unperturbed Hamihonian the result will 
depend on the choice of gauge in the one-photon ionisation transition-matrix (10). The 
length gauge in L^^^ gives similar results that LG SFA. On the other side, when L^^^ is 
calculated in the velocity gauge, the resulting spectra magnitudes are reduced to values 
similar to those obtained by solving numerically the TDSE. 

In figure 5 we compare the TDSE results with those for LG SFA, and with DipA- 
VG (SFA) results (8). For the DipA factorization, we employed plane- waves in the 
final-state and velocity gauge i.e. the right-side of (10). 

We observe that LG SFA spectra are at least one order of magnitude higher than the 
exact results while the velocity gauge in DipA-VG reduces considerably this difference, 
specially at the first peak where there is the major contribution to total ionisation 
probabilities. For the highest frequency investigated, DipA-VG gives almost the same 
spectra that TDSE while LG SFA remains one order of magnitude larger. 

There has been long standing discussions about the accuracy of different gauges in 
SFA [51, 52], that have recovered interest in recent years (see for example [53-56] to 
cite just a few), and deserves to be treated in a separated work. Nevertheless, here we 
simply note that in the multiphotonic regime, in the region of applicability of DipA, 
SFA calculations could be improved by choosing the velocity gauge in the computation 
of the factor L^^\ 
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Figure 6. (colour online) Ionisation spectra for H from Is solid line, 2s dotted line 
and 2po dashed dotted line for laser pulses with w = 0.855 a.u., Fq = 0.05 a.u., N = 7 
and 27 cycles. The ionisation spectra from 2s and 2po are shifted to lower energy in 
0.375 a.u. 



4.4- Ionisation from excited states and 2^^ -order approximation 

In this section we want to analyse the spectra for ionisation of hydrogen atoms from 
excited states. In figure 6 we present CVA ionisation probabihties for H atoms from 
the Is, 2s and 2po states. For the excited state we shifted the spectra in energy 
El — E2 = —0.375 a.u. As a result, the energy position of the ATI peaks for the 
different spectra are matched. We observe that the shape of the spectra are very similar 
in all cases, but the ionisation probabilities are lower for excited states. Theses facts 
can be understood in the DipA context: as we discussed before, the effect of the target 
on M^^^ produces a shift in the spectra that we have just corrected. Additionally, the 
intensity of the spectra is modulated by the different one-photon transition-matrices 
L^^^ in each case. 

Also, in figure 7 we show the ionisation probabilities from H(2s), comparing the 
CVA results with DipA. Like in the Is case (see figure 1), DipA underestimates the 
magnitude of ATI peaks for multiple photon absorption. Evidently, the first non- 
vanishing term in the multipolar expansion of the T-matrix is not enough to describe 
the ionisation process. In order to achieve a better representation we include a second 
term, i.e. T ?» T^^^ + T^'^' (see (7a)) in the differential ionisation probability (4). This 
calculation is also presented in figure 7, labelled as second-order (2'^'^T). We observe 
that the addition of the second term is sufficient to reproduce accurately the spectrum 
in an extended energy range. Visible differences appear only at high electron energies, 
at the third ATI peak. 
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Figure 7. (colour online) Ionisation spectra for H from 2s for laser pulses with w = 
0.855 a.u., Fo = 0.05 a.u. and N = 7 cycles. Full line (green): CVA, dotted line 
(black): DipA and dashed dotted line (blue): 2"^^ order calculations. 



5. Spectra for H2' 

In this section we study the ionisation of fixed-in-space H^ molecular ions by short 
laser-pulses. The nuclei of the molecule, having charges Z^ = Z;, = 1, are fixed at the 
internuclear distance R = 2 a.u. Employing spheroidal coordinates and using standard 
computational methods we calculate exactly the initial bound (pi{r) [57], and final 
continuum ^j{r) [58, 59] states of the molecule. Also the plane- wave exp {—iA{t) ■ r) is 
expressed in spheroidal coordinates [44, 60]. With these wave functions the factor g{t) 
in {6d) is evaluated numerically for the CVA calculation. To obtain DipA spectra, X^^^ 
is computed in the same way that in previous works [44, 61, 62]. 

In figure 8 we present the ionisation probability differential in energy and angle, for 
forward configuration, i.e: the emission direction is parallel to the internuclear axis and 
both are parallel to the polarisation vector (see figure 9(a)). As before, there is good 
agreement between the CVA and DipA calculations for all the energy-range investigated. 

Overall, the spectra for molecular ionisation are similar to the ones presented for 
atoms (figure 1 and 2), but there are some notable differences: 

(1) The spectra are energy-shifted due to different ionisation potentials of the targets. 
The computed Isa electronic energy is —1.1 a.u. -|- 1/R = —16.3 eV. Since we are 
considering fixed nuclei, the term 1/R is also present in the energy of the final 
continuum state. Then, the net shift in the spectra respect the atomic case is 
—0.6 a.u. as it can be seen in the position of ATI peaks in the figures. 

(2) The molecular spectra are shown for forward-direction emission, while in the atomic 
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Figure 8. (colour online) Electronic spectra for ionisation of HJ(Isct) as function of 
the electron energy, for emission in the forward direction. We have considered laser- 
pulses of A^ = 1,7, and 27 cycles, with frequencies uj = 0.855 and 1.71 a.u., and an 
amplitude Fq = 0.05 a.u. Full line: CVA and dotted line: DipA. 
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case we have presented the probabihties integrated in emission direction. The 
principal difference is that the deep minima in the forward condition, similar to the 
atomic spectra for fixed direction k \\ e (not shown), are softened by integration on 
emission angles. 

(3) The agreement between CVA and DipA calculations is better in the atomic case. 
This fact can be understood since the Isa orbital is more extended in space than 
the atomic Is. Then, for the same laser pulse, the approximation A ■ r <^ 1 is not 
satisfied in all range of the integral g{t). 

To analyse the emission in other directions we present in figure 9 the ionisation 
probabilities in three geometrical configurations: (a) in the forward direction, where 
the three vectors are parallel, i.e. Or = 9^ = 0, (b) with the electronic momentum 
parallel to the internuclear axis and both are at 45° respect to the polarisation vector, 
i.e. Or = Ok = 7r/4, and (c) where the molecule is perpendicular to i, Or = 90° and 
Ok = 45°. In all theses cases the spectra are well reproduced by DipA (not shown). In 
particular, the arrangements (b) and (c) have the same contribution oi i ■ k. Then, the 
factor M^^^ is exactly the same in both cases, and all dependence of the spectra with 
the molecular orientation arises from the one-photon ionisation transition-matrix. 

Similar conclusions can be deduced from the analysis of figure 10, where we compare 
the CVA results for ionisation from Isa, 2pa, in both cases with R = 2 a.u. and from Isa 
at -R = 1.4 a.u. This value is the equilibrium internuclear distance of the ground-state 
for the H2 molecule. As before, the differences are originated in the factor L^^\ 

As concluded from the analysis on molecular ionisation, the dependence on the 
emission direction, molecular orientation, and molecular structure are derived from the 
one-photon transition matrix properties. These results are consistent with reports on 
several works of laser-pulse ionisation of H2 and H^ where the authors have observed 
coincidences with the one-photon case [63-65]. Our present work shows that the 
possibilities of finding "new molecular effects" due to the interaction with laser pulses 
should be investigated outside the range where DipA reproduces the exact spectra. 

6. Conclusions 

We have investigated ionisation of atomic and molecular ions of hydrogen, induced 
by short laser pulses from ground and some excited states. Differential and total 
ionisation probabilities were calculated using the Coulomb- Volkov (CVA) and DipA 
approximations proposed in this work. The accuracy of the proposed approximations 
was studied by comparing with numerically exact solutions of the TDSE. The present 
results show that the calculations accurately describe the ionisation probabilities in the 
range of laser parameters here explored. 

DipA is deduced from the multipolar expansion of CVA T-matrix and therefore 
can be applied in the regime of validity of CVA, when is verified the condition u > Ip, 
while simultaneously being in the multiphotonic regime (F ^ 1). The main feature of 
DipA is that it expresses the ionisation spectra as a product of two factors. One of them 
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Figure 9. (colour online) CVA results for H^ ionisation spectra for the three 
geometrical arrangement: (a) forward configuration in solid line, (b) i? || fc at 45° 
of e in dotted line and (c) i? _L e and k at 45° in dashed dotted line. The laser 
parameters are the same than in the bottom panel of figure 8. 



accounts for the target structure and the ejected electron and is the well-known one- 
photon transition-matrix. The second factor depends comprises all the dependence on 
the laser-electron interaction and describes the main characteristics of the ATI spectra. 
Then, in this approximation laser and target contribute separately to the observed 
spectra. 

We have proposed several applications of DipA, covering several aspects of the laser- 
matter interaction, such as the interpretation of the features of the ATI spectra that 
come from each contribution. Also, the proposed approximations to first and second 
order are concrete examples of computationally low-cost alternative methods to calculate 
differential and total cross-sections in complex atomic and molecular targets. In all cases 
we discussed the limits of validity and range of applicability of the approximations. 
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Curves (b) and (c) are shifted in energy. 
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